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In recent work a general solution of the Ornstein Zernike equation for a general Yukawa closure 
for a single component fluid was found. Because of the complexity of the equations a simplifying 
assumption was made, namely that the main scaling matrix T had to be diagonal. While in principle 
this is mathematically correct, it is not physical because it will violate symmetry conditions when 
different Yukawas are assigned to different components. In this work we show that by using the 
symmetry conditions the off diagonal elements of T can be computed explicitly for the case of two 
Yukawas, and that although the solution is different than in the diagonal case, the excess entropy 
is formally the same as in the diagonal case. Analytical expressions for the Laplace transforms of 
the pair distribution functions are derived. 

PACS numbers: 61.20.Gy 



I. INTRODUCTION 

There are many problems of practical and academic 
interest that can be formulated as closures of some kind 
of either scalar or matrix Ornstcin-Zernike (OZ) equa- 
tion. These closures can always be expressed by a sum 
of exponentials, which do form a complete basis set if we 
allow for complex numbers Q, • 

While the initial motivation was to study simple ap- 
proximations like the Mean Spherical (MSA) or General- 
ized Mean Spherical Approximation (GMSA), the avail- 
ability of closed form solutions for the general closure 
of the hard core OZ equation makes it possible to write 
down analytical solutions for any given approximation 
that can be formulated by writing the direct correlation 
function c(r) outside the hard core as 



A I 



A I 



(n) f 



>/r = 



7r (1) 



In this equation is the interaction/closure constant 
used in the general solution first found by Blum and Hoye 
(which we will call BH78) 0, while K.W is the definition 
used in the later general solution by Blum, Vericat and 
Herrera (BVH92 in what follows) |g. In this work we 
will use the more common notation of BVH92. The case 
of factored interactions discussed by Blum, Q was sim- 
plified by Ginoza [lfjl f ll|. IT3 . IT^ | who found that as in 
the case of electrolytes [l4| the solution of the one expo- 
nent case could be expressed in terms of a single scaling 



parameter V . In the factorizable case we have 



where we have defined 



6^ = >) e - W2 



(3) 
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The general solution of this problem was formulated 
in by Blum, Vericat and Herrera |8( in terms of a scaling 
matrix T. The full solution was given recently by Blum 
et al. P, [Til Il6j| . For only one component the matrix 
r was assumed to be diagonal diagonal and explicit 
expressions for the closure relations for any arbitrary 
number of Yukawa exponents M were obtained. The 
solution is then remarkably simple in the MSA since 
then explicit formulas for the thermodynamic properties 
are obtained. 

The diagonal assumption is however not correct for 
mixtures, even if they are of the same hard core diameter. 
In this work we use the symmetry relations to calculate 
explicitly the off diagonal terms of T in the 1 component, 
2 yukawa case. 



II. SUMMARY OF PREVIOUS WORK 

We study the Ornstein-Zernike (OZ) equation 

hij(12) = c tJ (12) + J2 I d3h lk {lZ) Pk c k] {32) (4) 

k J 

where /iy(12) is the molecular total correlation function 
and Cij (12) is the molecular direct correlation function, 
Pi is the number density of the molecules i, and i = 1, 2 



2 



is the position f\ , r\i = \r\ — ri\ and cr.y is the dis- 
tance of closest approach of two particles i, j. The direct 
correlation function is 



M 



n=l 



r > cr,; 



(5) 



and the pair correlation function is 

M r ) = -! = -!> r < o-ij (6) 

We use the Baxter- Wertheim (BW) factorization of the 
OZ equation 



I + pH(k) I-pC(fc) =1 



(7) 



where / is the identity matrix, and we have used the 
notation 



HO) = 2 / drcos{kr)J(r) 
Jo 



C(fc) = 2 / drcos(fcr)S(r) 
Jo 

The matrices J and 5 1 have matrix elements 



(8) 



(9) 



Similarly, from Eq. I|12|) and Eq. Q we get, using the 
analytical properties of Q and Cauchy's theorem 

J(r) = Q(r) + / drxJ(r - r 1 )pQ(r 1 ) (16) 

The general solution is discussed in and yields 



A I 



q,., 



Xji < r (17) 



q%{r) = {l/2)A[{r-a/2f-(a l /2f] + 
p[(r - a/2) - fo/2)] + 

M 

£3 



n=l 



(n) e -z n a/2[ e -z n (r-a/2) _ e -z„<7i/2j. 

Xji < r < Uij 



(18) 



From here 



or 



jr i (n) - o - i 0o(^)ni n) = 



(19) 



dsshij(s) 



(10) 



2 * {-4" )n ^ + } = 5(n) (2°) 



Sij(r) — 2ir / dssCij(s) 



(11) 



I-pC(fc) = I-pQ(fc) I- i oQ T (-fc) (12) 



where Q T (— fc) is the complex conjugate and transpose 
of Q(fc). The first matrix is non-singular in the upper 
half complex fc-plane, while the second is non-singular in 
the lower half complex fc-plane. 

It can be shown that the factored correlation functions 
must be of the form 



Q(k) = I - p I dre lkr Q{r) 



where we used the following definition 



(13) 



(14) 



S(r) = Q(r) - J dr 1 Q(r 1 ) P Q T (r 1 - r) (15) 



A. The Laplace Transforms 

From Eq. we obtain the Laplace transform of the 
pair correlation function 

2n^2gn>(s)[5 ej - peqej(is)} = <?° (is) (21) 



where 



q%{is)= / dre~ sr [g?.(r)]' 



[(l + fp)A i+ ^ 

£- 

Z — / _L 



The Laplace transform of Eas. ltTTjl and l(TB|) yields 



(22) 



(23) 



£ — 

^ s + z, 



-[(Clf+Dlf)e-*^ 



i 2mO-^ (sCr i )C'i J - '' 
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This result will be used below. 

Another important relation deduced from Eg.|21jbv 
setting 



we write 



» 



y ' A nm A, 



(m) 



(33) 



An) 



(24) 



where 



Mr, 



(25) 

III. THE GENERAL CLOSURE 



The closure relation (BVH92 8] ) is, for only one com- 
ponent 

2TTKS {n) /z n + a^I^ - " pa {n) a {m) 

m 

j{n) p(m) _ Zm X( m )] - X^X^A = (26) 

For the one component case this simplifies to 
\^K6 nm + z n A (nm ^ + Zn [pa^a^A = 



where, as was shown in reference 8], A must be a 
symmetric matrix. 

From the symmetry of the contact pair correlation 
function Eq. I|16[l we get 

{ffy Oy ) = gji ((Jij } => { 9ij (vy)' = qjiivij)'} (34) 
which are 

^(n< n) - z„Af l V (n) = E( n(n) - ^ (n) k (n) (35) 



from which we get the scaling relation 

and a new set of M(M — l)/2 symmetry relations 



(36) 



(37) 



This can also be written as 



(27) 



27rA>[A (ll) ] 2 + z n pa (n) X {n) 
Y,^ 1 — [pa (n) a (m) ] [pA( m >AM = (28) 



Furthermore, using the scaling relations we get 

M A = T (38) 
where the matrix M (see Eq.j2Hl)has elements 

[M]„ m = — ^> \4 n ) {zmX ( T >} + Af^nf" 

(39) 

Solving these equations yields the relations 



which is the desired expression. This equation simplifies 
to 



2itpK n 



x (n) 



■ Z n (3 



{n) 



-/3 



(m) 



and 



where /3^ n ^ is 



= 
(29) 

(30) 



ivr 1 r = a 



(i + z-r -1 ) m = T 



(40) 



(41) 



IV. SYMMETRY 

In this section we will summarize and extend ourpre- 
vious analysis of the most general scaling relation |8j for 
the multiyukawa closure of the Ornstein Zernike equa- 
tion. We have 

ni w) = -X;r nTn x< m) 



(31) 



where T mn is the M x M matrix of scaling parameters. 
This matrix is not uniquely defined by the MSA closure 
relations and must be supplemented by M(M — 1) equa- 
tions obtained from symmetry requirements for the cor- 
relations. From the symmetry of the direct correlation 
function at the origin, Eq. (fTK|l 

qij(Xji) = qji(Xij) (32) 



Both T and A must be symmetric matrices. We have 
therefore a total of M(M — 1) symmetry relations, which 
together with the M closure equations give the required 
equations for the M 2 elements of the matrix T. 
The symmetry requirements are more explicitly 

r-M T = M-r T M T -[r T ] _1 = r _1 -M 5/(42) 

and 

(i + z • r^ 1 ) • m = m t • (1 + [r- 1 ] T • z) SII (43) 

the matrix M as 

M = - f) + ^M A (44) 

and 
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[M A U = — Sft k (n) (z m ^ m) - 2II< m >) - (z„^ (n) - 2Ilf ) )X<f l) ' 

^— Y> fx(")xW( 2m <5 pm + 2r mp ) - (^5 pn + 2r np )x( m )xW 



[M A ] nm = -X^X^ [ lnm + , 



(45) 



where 



2p 



■E 



r r 



and 



7* 



The second symmetry condition is Ea. (|43|) is 
M A = T 1 ■ M ■ z - z ■ M T • [r T ]- x 



(47) 



V. THE 2 YUKAWA CASE: SYMMETRIC 
MATRIX RESULTS 

We write equation l|24|l in matrix form |16| 

-5, =M-a, 

where 



(i) 

i 

(2) 



n, 



(i) 

i 

(2) 



4 15 

,( 2 ) 



Using the symmetry relation ea. (|42|l we get 



(50) 



with 



X12 = 



712 



21 - Z2 



zi + z 2 + 2r( n ) +2r( 22 ) 
z 1 -z 2 + 2r( 11 )-2r( 22 ) 



Zl + 2 2 



(52) 



(53) 



in ea. H46(l we can write 



I 

(46) We rewrite eq.lgSJ as 
2 



- 0(1) - 




. s (2) . 





1 1 - Xl2 
1+X12 1 



/3W 



(55) 



Here we have defined 

nil) - r( n ) + ^.r( 12 >- 
y XW 

If we also define 

(48) 

then we can solve eq.(JSSJ) 
20« =2/3 s + /? 12X i 2 ; 



ft = -2 



X12 



C (2) _ p(22) , ^ll r (21) 

(56) 

2 g(i2) = g0) _ g(a) ( 57 ) 



2<? (12) = -AX12; (58) 
g(i2)- 



/3i2 = 



2 

X12 



X12 



(59) 



(49) From the second symmetry condition Ea. (|60|l we get 

Vi2) Zl r(i2) _ ^ Z2 T^ - 2T^T^ X i2 + 2r 12 = 

(60) 

(61) 



where 

TL2 = 



X(2)' 



z 2 r( n ) { Zl + r^ 11 )) - Zl r(22) ( Z2 + r( 22 )) \ 



zi + z 2 + 2r( n ) + 2r( 22 ) 

We remark also that 

7i2 = \[z2T (11) (i + xi 2 ) - zir( 22 )(i - X12)] 



) 



+«2r< 111 r< 22 > (62) 



in Ea. (|6T)j> we get 

fS^ r(12) -fS^) + 2 Xl2 ^ r 

+[z 2 r( 11 )(i + Xl2 ) - Zi tW(i - X12)] = (63) 

Using now ea. (|5*T)l 



X« 

x^ 



1 1-X12 
1 + X12 1 



xw 
x^ 



X 2 r 2 

x\ 12 



X2- 

Xx 



S12X12 



Z2 



2 p(22) 



(54) 



+ r( 22 )(z 2 + r( 22 >) = 

(64) 



5 



and 



xf 21 ~ x? 21 



JJ2X12 +zi + 2r (ll) 

+r( 11 )(z 1 + r( 11 )) 



o 



from where 



e '" - 5 [{ 



S12X12 . Zl2 I / A (l) 



with 



A r } = Z 2 + s 12Xl2 



S12X12 



(22) 



a[! 5 = z\ + S12X12 
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S12X12 



+ z 2 + 2r( ; 
- zi - 2r< n > 



(65) 
(66) 
(67) 

(68) 
(69) 



2 ^-{f + - } 



X12 + z\ 



2g m _ {_£!?. + Zl2 } Xl2 + Z2 

From here we get the equation 



2^12 v J 

X12 - T~nrT } i X12 



z 12 + 2g^ 



S12 + 2Q S J L "1.2 
which yields the two solutions 

2^12 



(70) 
(71) 



X12 



S12 + 2Q S 



(A); xi 



£12 + 2g< 12 ) 
S12 



(72) 



Notice first that in the zero density limit we get 

X12 =>■ ; X12 — (73) 

S12 S12 

and then in eas H66(l and (|67|l we get the correct zero den- 
sity limit only from choice (B) 



gW ~ L 
2 



2 



( Z12 . \ S12 

{— + 'uj-**- T 



(74) 
= (75) 



Then 



012 = 2- 



8 --2G {12) — ■ 



S12 



Z12+2GW) 



lis) 



2G^s 12 



( Zl2 + 2gW) 



(76) 



These expressions turn out to be identical to those de- 
rived by Blum and Ubriaco using the diagonal approxi- 
mation Q] 



VI. THERMODYNAMICS BY PARAMETER 
INTEGRATION 

We will use the notation and results of Blum and Her- 
nando 111 . We recall that 



j-(n)jj(n) _ j(n)j£"(n) _ fi(n) 

Remember that 
Here 

= 6irMzn°) ~ 2pl3°a 3 Mzn) 

and 

z n g 

The total excess internal energy is 
E(/3) 



kTV 

n 

From ea. (|31|l we show that 

_tt(") = gWjW 



(77) 
(78) 
(79) 

(80) 

(81) 
(82) 



where Q^ n > is a (generally algebraic) function of the co- 
efficients f3 = /9a 1 ■■}■ I n fact m eq.JTTJ) 



= ^{J (n) <5 nm + j(«) r («™)}x(" 1 ) 

m 

. j(n)j£-(n) _|_ j-(n) p(nm)|j^(m) 
m 



with 



For the 1 component case we get 



J(n) _|_ g(n)j(n) 



(83) 



(84) 



(85) 



Then, since the 'charge' parameters are constants at 
constant temperature, the derivative of flW with respect 
to the scaling parameter is 



0jj(n) 



J 



(n) 



J 



(n) 



f d (IW) 
I <%?(«) 



(J(n) + g(n) 



(86) 
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where we use the fact that J^ n > is independent of Q^ n \ 
The desired energy derivative Eg. 1)8 1|) are 



dE 



= -p[X^} 2 



or 



dE 



dE 



= - P {[X^} 2 - [X^] 2 } 



Qg(nrn) 

The integrability condition is satisfied since 
d 2 E d 2 E 



= 6: 



K 



j(n) 



j(n) _|_ g(n)j{n) 



(87) 



We use now Ea. (|30[l to obtain 

dE 1 r 9 , 
= - [P 2 S + s 12 /3 s + * 12 /? ia ] 



2(20(i 2 ) +z 12 ) 



20 (s) + s 12 



S\ 2 Z\ 2 



(20(12) +2l2 ) 



(90) 



(91) 



and 



dE 1 



dgw 2 

S12Z12 



P*(Pa +S 12 ) + ZviPs + ^{Pl ~ Pl 2 } 
2si2 



4(2gd 2 ) +2l2 )2 



s 2 12 + zf 2 -{2g^ + s 12 2si2Zl2 



(20d2) + ZJ3 

S12Z12 



(92 



(89) Thermodynamic integration of these equations leads to 



S12Z12 



(z 12 + 20d2)) 




1 - 



(z 12 + 20(i2))( Sl2 + 20( s ); 



S12Z12 



S12Z12 



(z X2 + 20d2)) 



AS 1 

-z 12 + 20< 12 ) 



.?? 



2tt 



12 



(93) 



AS = -- 



2vr 



Although the derivation of this equation is completely 



different from that obtained using the diagonal T 
assumption, the resulting entropy is identical to this 
case when proper reference states are used 
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